Interface-Roughening Phase Diagram of the Three-Dimensional Ising Model for 
All Interaction Anisotropies from Hard-Spin Mean-Field Theory 
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The roughening phase diagram of the d = 3 Ising model with uniaxially anisotropic interactions 
is calculated for the entire range of anisotropy, from decoupled planes to the isotropic model to 
the solid-on-solid model, using hard-spin mean-field theory. The phase diagram contains the line 
of ordering phase transitions and, at lower temperatures, the line of roughening phase transitions, 
where the interface between ordered domains roughens. Upon increasing the anisotropy, roughening 
transition temperatures settle after the isotropic case, whereas the ordering transition temperature 
increases to infinity. The calculation is repeated for the d — 2 Ising model for the full range of 
anisotropy, yielding no roughening transition. 
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I. INTRODUCTION 

The ordering phase transition in a crystal precipitates 
the formation of macroscopic domains, differently or- 
dered with respect to each other. The interface between 
such domains incorporates static and dynamic phenom- 
ena of fundamental and applied importance. Of singu- 
lar importance is the occurrence of yet another phase 
transition, distinct from the ordering phase transition, 
which is the interface roughening phase transition. [l|, 0| 
The roughening phase transition is well-studied with the 
three-dimensional Ising model, in the so-called solid-on- 
solid limit, in which the interactions along one spatial 
direction (z) are taken to infinite strength, while the in- 
teractions along the x and y spatial directions remain 
finite. In this case, due to the infinite interactions, the 
ordering phase transition moves to infinite temperature 
and is not observed. A study of the system with finite 
interactions, where both ordering and roughening phase 
transitions should distinctly be observed, had not been 
done. 

In our current study, hard-spin mean- field theory 0,0], 
which has been qualitatively and quantitatively success- 
ful in frustrated and unfrustrated magnetic ordering 
problems pj-flU , is used to study ordering and roughening 
phase transitions in the three-dimensional (d — 3) Ising 
model for the entire range of interaction anisotropies, 
continuously from the solid-on-solid limit to the isotropic 
system to the weakly-coupled-planes limit. The phase 
diagram is obtained in the temperature and interaction 
anisotropy variables, with separate curves of ordering and 
roughening phase boundaries. The method, when ap- 
plied to the anisotropic d = 2 Ising model, yields the 
lack of roughening phase transition. 



II. HARD-SPIN MEAN-FIELD THEORY 

Hard-spin mean-field theory has been introduced as a 
self-consistent theory that conserves the hard-spin (|s| = 
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FIG. 1: (Color online) For the d = 3 anisotropic Ising model, 
magnetizations mi versus xy layer-number i curves for differ- 
ent temperatures 1/J xy . Each panel shows results for the in- 
dicated anisotropy J z /J xy . The curves in each panel, with de- 
creasing sharpness, are for temperatures 1/J xy ~ 1, 3, 5, 6. In 
the upper panels, the high-temperature curves coincide with 
the horizontal line mi — 0. 



1) condition, indispensable to the study of frustrated 
systems. 0, 0] This method is almost as simply imple- 
mented as usual mean-field theory, but brings consider- 
able qualitative and quantitative improvements. Hard- 
spin mean-field theory has yielded, for example, the lack 
of order in the undiluted zero-field triangular-lattice an- 
tifcrromagnetic Ising model and the ordering that occurs 
either when a uniform magnetic field is applied to the 
system, giving a quantitatively accurate phase diagram 
in the temperature versus magnetic field variables 0- 
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FIG. 2: (Color online) Local magnetization data for the 
d = 3 anisotropic Ising model. The curves, starting from 
the high-temperature side, are for anisotropics J z /J xy = 
10, 5, 2, 1, 0.5, 0.2 . Upper panel: Magnetization absolute val- 
ues \m.b\ away from the interface as a function of temperature 
1/Jxy, for different values of the anisotropy J z /J xy . Lower 
panel: The deviation |mj,| — |m, | averaged over the system ver- 
sus temperature 1/J xy for different anisotropies J z /J xy . This 
averaged deviation vanishes when the interface is smooth. 
Note the qualitatively different low-temperature behavior in 
the d = 2 case shown in Fig. 4. 



@, S o r when the system is sublattice-wise quench- 
diluted [17|- Hard-spin mean-field theory has also been 
successfully applied to complicated systems that exhibit a 
variety of ordering behaviors, such as three-dimensionally 
stacked frustrated systems d, 0] , higher-spin systems Q , 
and hysteretic d = 3 spin-glasses [l8j]. Furthermore, 
hard-spin mean-field theory shows qualitative and quan- 
titative effectiveness for unfrustrated systems as well, 
such as being dimensionally discriminating by yielding 
the no-transition of d — 1 and improved transition tem- 
peratures in d — 2 and d = 3.0, [l8j 

We have therefore applied hard-spin mean-field theory 
to the global study of the roughening transition in the 
anisotropic d = 3 Ising model. [We have also found that 
no roughening phase transition is seen in d = 2 (Sec. IV)]. 
The uniaxially anisotropic d = 3 Ising model is denned 
by the Hamiltonian 



xy 



'li = Jxy > *;* 



Jz 



(1) 



where, at each site i of a d = 3 cubic lattice with pe- 
riodic boundary conditions, Si = ±1. The first sum is 
over nearest-neighbor pairs of sites along the x and y 



spatial directions and the second sum is over the nearest- 
neighbor pairs of sites along the z spatial direction. The 
interactions are ferromagnetic, J xy ,Jz > 0, except for 
the interaction between two of the xy planes, which has 
the same magnitude as the other J z interactions, but 
is antiferromagnetic: J z = — J z < 0. This choice is 
made in order to induce an interface when the system 
is ordered. (An alternate approach would have been to 
use a system without periodic boundary conditions along 
the z direction, but with oppositely pinned spins at each 
edge. However, this would have introduced a surface ef- 
fect at the pinned edges, modifying the magnetization 
deviations which would thereby not exclusively reflect 
the spreading of the interface.) 

For this system, the self-consistent equation of hard- 
spin mean-field theory is 



in,. 



E 




x tanh 



where the last sum is over the sites j that are nearest 
neighbor to site i and the first sum is over all states {sj} 
of the spins at these nearest-neighbor sites. In Eq.(2), 



P{m j:Sj ) = - (1 + mjSj) 



(3) 



is, for local magnetization rrij at site j, the probabil- 
ity of having the spin value of Sj. The coupled Eqs. 
(2) are solved numerically for a 20x20x20 cubic system 
with periodic boundary conditions, by iteration: A set of 
magnetizations is substituted into the right-hand side of 
Eqs. (2), to obtain a new set of magnetizations from the 
left-hand side. This new set is then substituted into the 
right-hand side, and this procedure is carried out repeat- 
edly, converging to stable values of the magnetizations 
that is the solution of the equations. The resulting mag- 
netization values depend on the z coordinate only. 



III. RESULTS: ORDERING AND 
ROUGHENING PHASE TRANSITIONS IN d = 3 

A series of curves for the magnetizations mi versus 
xy layer number i are shown for different temperatures 
1 / J xy , for a given anisotropy J z / J xy in each panel of Fig. 
1 . For each value of the anisotropy, the magnetizations 
rrii are zero at high temperatures and become non-zero 
below the ordering transition temperature Tq- The or- 
dering onset is seen in the upper panel of Fig. 2, where 
the magnetization absolute values \rti\\ away from the in- 
terface are plotted as a function of temperature l/J X y, 
for different values of the anisotropy J z / J xy . 

In Fig. 1, it is also seen that, at temperatures just 
below Tc, the interface between the rrij ^ domains 
is spread over several layers. It is also seen that be- 
low a lower, roughening-transition temperature Tr, the 
interface becomes localized between two consecutive lay- 
ers, reversing the sign of the magnetization with no 
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FIG. 3: (Color online) For the d — 3 anisotropic Ising model, 
the calculated phase diagram showing the disordered, or- 
dered with rough interface, and ordered with smooth interface 
phases. The squares indicate the exact ordering temperatures 
from duality at J z / J xy — and from Ref. at J z / J xy = 1. 
The circle indicates the roughening transition temperature 
for the solid-on-solid limit J z /J xy — > oo.Q] The roughening 
transition is obtained by fitting the averaged deviation curves 
(lower panel of Fig. 2) within the range < |m;,| — \m,i\ >= 0.01 
to 0.04, to find the temperature at which the averaged devi- 
ation reaches zero, meaning that the interface becomes local- 
ized between two consecutive layers, reversing the sign of the 
magnetization m b with no change in magnitude. 



Thus, we have deduced the phase diagram, for all val- 
ues of the anisotropy J z /J xy and temperature 1/J xy , as 
shown in Fig. 3. The roughening transition is obtained 
by fitting the averaged deviation curves (lower panel of 
Fig. 2) within the range < \m b \ - \rrn\ >= 0.01 to 0.04, 
to find the temperature at which the averaged deviation 
reaches zero, meaning that the interface becomes local- 
ized between two consecutive layers, reversing the sign 
of the magnetization m,b with no change in magnitude. 
In Fig. 3, the ordering and roughening phase transitions 
occur as two separate curves, starting in the decoupled 
planes (J z / Jxy = 0) limit and scanning at finite temper- 
ature the entire range of anisotropies. The ordering tran- 
sition starts, for the decoupled planes limit J z /J xy = 0, 
at 1/J xy =3.12, to be compared with the exact result 
of l/J X y — 2.27. The ordering transition continues to 
1/Jxy = 5.06, to be compared with the precise [lj| re- 
sult of l/J xy = 4.51, for the isotropic case J z /J xy = 1. 
In the solid-on-solid limit (J z /J xy — ► oo), the ordering 
boundary goes to infinite temperature. The roughening 
transition starts at 1/ J xy = for J z / J xy close to zero and 
settles to a finite temperature value before the isotropic 
case. Thus, the roughening transition temperature 1/ J xy 
is 1.45 in the isotropic case J z /J xy = 1 and 1.62 in the 
solid-on-solid limit J z /J xy — > oo, the latter to be com- 
pared with the value of 2.30 ± 0.10 from Ref.0. 



IV. RESULTS: ORDERING TRANSITIONS BUT 
NO ROUGHENING TRANSITIONS IN d = 2 



change in magnitude. This onset is best seen in the lower 
panel of Fig. 2, where the deviation \nib \ — averaged 
over the system is plotted as a function of temperature 
1/ 'Jxy for different anisotropies J z /J xy . 
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FIG. 4: (Color online) For the d — 2 anisotropic Ising model, 
the deviation \m^\ — \ rrii \ averaged over the system versus tem- 
perature l/Jxy for different anisotropies J z /J xy . The curves, 
starting from the high-temperature side, are for anisotropies 
Jz/Jxy = 10, 5, 2, 1, 0.5, 0.2 . It is seen that the deviation does 
not vanish, i.e., the interface does not localize, down to zero 
temperature. Thus, a qualitatively different low-temperature 
behavior occurs, as compared with the d — 3 case shown in 
the lower panel of Fig. 2. 



We have also applied our method to the anisotropic 
d=2 Ising model, defined by the Hamiltonian 



/3H = J x SiSj + J z 5 ' Sis 



(ij) 



(4) 



where, on a 20x20 square lattice with periodic boundary 
conditions, the first sum is over nearest-neighbor pairs 
of sites along the x spatial direction and the second sum 
is over the nearest-neighbor pairs of sites along the only 
other (z) spatial direction. 

The ordering phase transition is observed in d = 2 
similarly to the d = 3 case. However, the rough inter- 
face phase continues to zero temperature, as seen in the 
\mi,\ — \m,i\ curves in Fig. 4. Thus, no roughening phase 
transition occurs in d — 2. The corresponding phase di- 
agram is given in Fig. 5. The ordering transition starts, 
for the decoupled lines limit J z /J x = 0, at l/J x = 0, 
as expected for decoupled d — 1 systems. The order- 
ing transition continues to 1/J X =3.09, to be compared 
with the exact result of l/J x — 2.27, for the isotropic 
case J z /J x = 1. In the J z /J x — > oo limit, the ordering 
boundary goes to infinite temperature. 
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sitions for the isotropic and anisotropic Ising models, in 
spatial dimensions d = 3 and 2. This result attests to 
the microscopic efficacy of the model. Future works, 
such as the effects of uncorrelated and correlated (aerogel 
[20l [2l|) frozen impurities on the roughening transitions, 
are planned. 



FIG. 5: (Color online) For the d=2 anisotropic Ising model, 
the phase diagram showing the disordered phase and the or- 
dered phase with rough interface. The dashed curve is the ex- 
act ordering boundary sinh(2J x )sinh(2J z ) = 1 obtained from 

duality. No ordered phase with smooth interface is found. Acknowledgments 



V. CONCLUSION 

It seen that hard-spin mean-field theory yields a com- 
plete picture of the ordering and roughening phase tran- 
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